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Abstract 

This paper presents an analysis of target localization accuracy, attainable by the use of MIMO (Multiple-Input 
Multiple-Output) radar systems, configured with multiple transmit and receive sensors, widely distributed over a given 
area. The Cramer-Rao lower bound (CRLB) for target localization accuracy is developed for both coherent and non- 
coherent processing. Coherent processing requires a common phase reference for all transmit and receive sensors. 
The CRLB is shown to be inversely proportional to the signal effective bandwidth in the non-coherent case, but is 
approximately inversely proportional to the carrier frequency in the coherent case. We further prove that optimization 
over the sensors' positions lowers the CRLB by a factor equal to the product of the number of transmitting and 
receiving sensors. The best linear unbiased estimator (BLUE) is derived for the MIMO target localization problem. 
The blue's utility is in providing a closed form localization estimate that facilitates the analysis of the relations 
between sensors locations, target location, and localization accuracy. Geometric dilution of precision (GDOP) contours 
are used to map the relative performance accuracy for a given layout of radars over a given geographic area. 

Index Terms 

MIMO radar, spatial processing, adaptive array. 

I. Introduction 

A. Background and Motivation 

Research in MIMO radar has been growing as evidenced by an increasing body of literature [1]-[18]. Generally 
speaking, MIMO radar systems employ multiple antennas to transmit multiple waveforms and engage in joint 
processing of the received echoes from the target. Two main MIMO radar architectures have evolved: with colocated 
antennas and with distributed antennas. MIMO radar with colocated antennas makes use of waveform diversity [2], 
[4], [12], [14], [15], while MIMO radar with distributed antenna takes advantage of the spatial diversity supported by 
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the system configuration [1], [3], [5], [13]. MIMO radar systems have been shown to offer considerable advantages 
over traditional radars in various aspects of radar operation such as the detection of slow moving targets [17], 
[8], the abiUty to identify and separate multiple targets [10], [11], and in the estimation of target parameters such 
as direction-of-arrival (DOA) [8], [10], and range-based target localization [18]. In particular, [18] studies target 
localization with MIMO radar systems utilizing sensors distributed over a wide area. 

Conventional localization techniques include time-of-arrival (TOA), time-difference-of-arrival (TDOA), and direction- 
of-arrival (DOA) based schemes. MIMO radar system with colocated antenna can perform DOA estimation of targets 
in the far-field, in which case, the received signal has a planar wavefront. In this class of systems, extensive research 
has focused on waveform optimization. In [7], [14], [15] the signal vector transmitted by a MIMO radar system is 
designed to minimize the cross-correlation of the signals bounced from various targets to improve the parameter 
estimation accuracy in multiple target schemes. Some of the waveform optimization techniques suggested in [16] 
are based on the Cramer-Rao lower bound (CRLB) matrix [19], [20]. The CRLB is known to provide a tight 
bound on parameter estimation for high signal-to-noise ratio (SNR). Several design criteria are considered, such as 
minimizing the trace, determinant, and the largest eigenvalue of the CRLB matrix, concluding that minimizing the 
trace of the CRLB gives a good overall performance in terms of lowering the CRLB. In [9], a CRLB evaluation of 
the achievable angular accuracy is derived for linear arrays with orthogonal signals. The use of orthogonal signals 
is shown to provide better accuracy than correlated signals. For low SNR scenarios, the Barankin bound is derived 
in [10], demonstrating that the use of orthogonal signals results in a lower SNR threshold for transitioning into the 
region of higher estimation error. 

MIMO radar systems with widely spread antennas take advantage of the geographical spread of the deployed 
sensors. The multiple propagation paths, created by the transmitted waveforms and echoes from scatterers in their 
paths support target localization through either direct or indirect multUateration. With direct multUateration, the 
observations collected by the sensors are jointly processed to produce the localization estimate. With indirect 
multilateration, the TOAs are estimated first, and the localization is subsequently estimated from the TOAs. The 
observations and processing can also be classified as either non-coherent or coherent. The distinction between the 
two modes relies on the need for mere time synchronization between the transmitting and receiving radars in the 
non-coherent case, versus the need for both time and phase synchronization in the coherent case. Note that our 
coherent/non-coherent terminology is limited to the processing for localization. Thus, a transmitted signal may have 
in-phase and quadrature components, yet the localization processing is non-coherent if it utilizes only information 
in the signal envelope. In the sequel, we evaluate the performance of localization utilizing both coherent and 
non-coherent processing. 

MIMO radar systems belongs to the class of active localization systems, where the signal usually travels a round 
trip, i.e. the signal transmitted by one sensor in a radar system is reflected by the target and measured by the 
same or a different sensor. Traditional single-antenna radar systems, performing active range-based measurements, 
are well known in literature [21]-[25]. The target range is computed from the time it takes for the transmitted 
signal to get to the target plus the travelling time of the reflected signal back to the sensor. The range estimation 
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accuracy is directly proportional to the mean squared error (MSE) of the time delay estimation and is shown to be 
inversely proportional to the signal effective bandwidth [21]. A first study of the localization accuracy capabiUty 
of MIMO radar systems is provided in [18], where the Fisher information matrix (FIM) is derived for the case of 
orthogonal signals with coherent processing and widely separated antennas. The CRLB is analyzed numerically, 
pointing out the dependency of the accuracy on the signal carrier frequency in the coherent case, and its reliance on 
the relative locations of the target and sensors. In [18], it is observed that the CRLB is a function of the number of 
transmitting and receiving sensors, however an analytical relation is not developed. The high accuracy capability of 
coherent processing is illustrated by the use of the ambiguity function (AF). Active range-based target locaUzation 
techniques are also used in multistatic radar systems, proposed in [26]. The TOA of a signal transmitted by a single 
transmit radar, reflected by the target and received at multiple receive anteimas is used in the localization process. 
It is observed that increasing the number of sensors improves locaUzation performance, yet an exact relation is not 
specified. This paper addresses deficiencies in the literature by obtaining closed-form expressions of the CRLB for 
both coherent and non-coherent cases. 

Geolocation techniques has been the subject of extensive research. Geolocation belongs to the class of passive 
localization systems, where the signal travels one-way. Since these passive measurement systems employ multiple 
sensors, further evaluation of existing results for geolocation systems might provide insightful for the active case. 
In wireless communication, passive measurements are used by multiple base stations for localization of a radiating 
mobile phone. The localization accuracy performance is evaluated in [27]. It is shown that the localization accuracy 
is inversely proportional to the signal effective bandwidth as it does in the active localization case. Moreover, 
the accuracy estimation is shown to be dependent on the sensors/base stations locations. In navigation systems, the 
target makes use of time synchronized transmission from multiple Global Positioning Systems (GPS) to establish its 
location. In [28], [29], the relation between the transmitting sensors location and the target localization performance 
is analyzed. GDOP plots are used to demonstrate the dependency of the attainable accuracy on the location of the 
GPS systems with respect to the target. In an optimal setting of the GPS systems relative to the target position, the 
best achievable accuracy is shown to be inversely proportional to the square root of the number of parlicipating 
GPS. In the sequel, we apply the GDOP metric to evaluate the locaUzation performance of MIMO radar. 

B. Main Contributions 

The main contributions of this paper are: 

1) The CRLB of the target localization estimation error is developed for the general case of MIMO radar 
with multiple waveforms transmission. The analytical expressions of the CRLB are derived for the case 
of orthogonal waveforms with non-coherent and coherent observations. The non-coherent case is used as 
benchmark for evaluating the performance of the system with coherent observations. 

2) It is shown that the CRLB expressions for both the non-coherent and coherent cases can be factored into two 
terms: a term incorporating the effect of bandwidth and SNR, and another term accounting for the effect of 
sensor placement. 



September 24, 2008 



DRAFT 



4 



3) The CRLB of the standard deviation of the localization estimate with non-coherent observations is shown 
to be inversely proportional to the signals averaged effective bandwidth. Dramatically higher accuracy can 
be obtained from processing coherent observations. In this case, the CRLB is inversely proportional to the 
carrier frequency. This gain is due to the exploitation of phase information, and is referred to as coherency 
gain. 

4) Formulating a convex optimization problem, it is shown that synmietric deployment of transmitting and 
receiving sensors around a target is optimal with respect to minimizing the CRLB. The closed form solution 
of the optimization problem also reveals that optimally placed M transmitters and N receivers reduce the 
CRLB on the variance of the estimate by a factor MN/2. This is referred to as the MIMO radar gain. 

5) A closed form solution is developed for the BLUE of target localization for coherent MIMO radars. It 
provides a closed form solution and a comprehensive evaluation of the performance of the estimator's MSE. 
This estimator provides insight into the relation between sensors locations, target location, and localization 
accuracy through the use of the GDOP metric. Contour maps of the GDOP, presented in this paper, provide a 
clear understanding of the mutual relation between a given deployment of sensors and the achievable accuracy 
at various target locations. 

The rest of the paper is organized as follows: The system model is introduced in Section II. In Section III, the 
CRLB is derived for the general case of multiple transmitted waveforms. Analytical expressions are obtained for the 
cases of non-coherent and coherent observations with orthogonal signals. Optimization of the CRLB as a function 
of sensor location is provided in Section IV. The performance of two localization estimators is evaluated in Section 
V. To establish a better understanding of the relations between the radar geographical spread and the target location, 
the GDOP metric is introduced in Section V-D. Finally, Section VI concludes the paper. 

A comment on notation: vectors are denoted by lower-case bold, while matrices use upper-case bold letters. The 
superscripts "T" and "H" denote the transpose and Hermitian operators, respectively. Complex conjugate is denoted 
0*. Points in the x-y plane are denoted in upper-case X = {x, y) . 

II. System Model 

We consider a widely distributed MIMO radar system with M transmitting radars and N receiving radars. 
The receiving radars may be colocated with the transmitting ones or individually positioned. The transmitting 
and receiving radars are located in a two dimensional plane {x,y). The M transmitters are arbitrarily located at 
coordinates = {xtu, Utk) , k = 1, . . . , M, and the N receivers are similarly arbitrarily located at coordinates Rg = 
{xre, yri) , ^ = ■ ■ ■ ,N. The set of transmitted waveforms in lowpass equivalent form is Sk{t) , k = 1, . . . ,M, 
where J^\sk{t)\^ dt = 1, and T is the connmon duration of all transmitted waveforms. The power of the transmitted 
waveforms is normalized such that the aggregate power transmitted by the sensors is constant, irrespective of the 
number of transmit sensors. To simplify the notation, the signal power term is embedded in the noise variance term 
such that the SNR at the transmitter, denoted SNRj and defined as the transmitted power by a sensor divided by the 
noise power at a receiving sensor, is set that a desired level. Let all transmitted waveforms be narrowband signals 
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with individual effective bandwidth (Sk defined as Pi = ( /^^ p \Sk (/)|^ df^ / ( /^^ \Sk (/)|^ d/) , where the 
integration is over the range of frequencies with non-zero signal content Wk [21]. We further define the signals 
averaged effective bandwidth or rms bandwidth as P^ = jjJ2h=if^k ^^'^ normaUzed bandwidth terms as 
PRk = Pk/P- The signals are narrowband in the sense that for a carrier frequency of /c, the narrowband signal 
assumption implies Pi/ <C 1 and /3^/ <C 1. 

The target model developed here generalizes the model in [21J to a near-field scenario and distributed sensors. In 
Skolnik's model [21], the returns of individual point scatterers have fixed amplitude and phase, and are independent 
of angle. For a moving target, the composite return fluctuates in amplitude and phase due to the relative motion of the 
scatterers. When the motion is slow, and the composite target return is assumed to be constant over the observation 
time, the target conforms to the classical Swerling case I model. We now proceed to generalize this model to a 
target observed by a MIMO radar with distributed sensors. Assume an extended target, composed of a collection 
of Q individual point scatterers located at coordinates Xq = {xq, yq) , q = 1, . . . ,Q. The ampUtudes (q of the 
point scatterers are assumed to be mutually independent. The pathloss and phase of a signal reflected by a scatterer, 
when measured with respect to a transmitted signal Sk {t) , are functions of the path transmitter-scatterer-receiver. 
Let Tik (Xq) denote the propagation time from transmitter k, to scatterer q, to receiver £, 

nk {Xq) = ^ (^^J {Xtk - Xqf + {ytk - Vqf + \J {XrC ^ Xqf + {Vri - Vqij , (1) 

where c is the speed of light. Our signal model assumes that the sensors are located such that variations in the 
signal strength due to different target to sensor distances can be neglected, i.e., the model accounts for the effect 
of the sensors/target localizations only through time delays (or phase shifts) of the signals. The conmion path loss 
term is embedded in The baseband representation for the signal received at sensor i is: 

M Q 

(0 = XI ^1 {-J'^^fcnk {Xq)) Sk {t - Tik {Xq)) + Wt{t), (2) 

fe=l g=l 

where the term 2'jTfcTek {Xq) is the phase of a signal transmitted by sensor k, reflected by scatterer q located at Xq, 
and received by sensor £. Phases are measured relative to a conamon phase reference assumed to be available at the 
transmitters and receivers. The term we {t) is circularly symmetric, zero-mean, complex Gaussian noise, spatiaUy 
and temporally white with autocorrelation function (r). The noise term is set = l/SNRj, where SNR* is 
measured at the transmitter. SNRj is normalized such that the aggregate transmitted power is independent of the 
number of transmitting sensors. The SNR at the receiver, due to a scatterer with amplitude Cg, is SNR^ = |Cg|^SNRt. 
Signals reflected from the target combine at each of the receive antennas. For example, the resultant signal at receive 

anterma i is given by 

Q 

X CijSfe {t - Tik {Xq)) exp {-j2nfcnk {Xq)) « ('sk {t - m {X')) exp {-j2TTfcnk {X')) , (3) 

9=1 

where and {2-KfcTik {X')) are respectively the amplitude and phase given by 

1/2 




Cq COS {2Tr fcTek {Xq)) + V Cg sin {2Tr fcTtk {Xq)) 



(4) 
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and 

„ , , -1 E?=iC,sin(27r/em(X,)) 

27rjcT|fe (a ) = tan — . (5) 

E?=iC,cos(27r/em (X,)) 

In obtaining (3), we invoked the narrowband assumption Sk {t — rek {Xq)) ~ Sk{t — r^k {X')), for all scatterers, 
namely that the change in the lowpass equivalent signals across the target is negligible. It follows from this discussion 
that the extended target is represented by a point scatterer of amplitude C' and time delays r^k {X') , where all the 
quantities are unknown. 

While this target model is completely adequate for our needs, it is possible to extend it sUghtly, at httle cost. 
Assume a constant time offset error Ar at the receivers. Further, assume that the error is small such that it does not 
impact the signal envelope, but it does impact the phase. Then we can write the time delays r^k {X') = r^k (X) + Ar 
for some location X = {x, y) . The target model (3) can now be expressed 

C'sfe {t - Tik {X')) exp (-i27r/em {X')) « C^fc (i - r^fe {X)) exp {-j27rf,rek {X)) , (6) 

where ( = (^'e~^'^'^f''^'^ and the narrowband assumption was invoked once more. The composite target of (3) is 
then equivalent to a point scatterer of complex amplitude ( and time delays Tik {X) . For simplicity, the following 
notation is used: r^k = T^k {X). The signal model (2) becomes 

M 

n (i) = X] ^ (-i27r/cm) Sk {t -Tek) + we{t). (7) 

fc=i 

We define the vector of received signals as r = [ri,r2, r^]^ for later use. The radar system's goal is to estimate 
the target location X = {x, y) . The target location can be estimated directly, for example by formulating the 
maximum likelihood estimate (MLE) associated with (7). Alternatively, an indirect method is to estimate first the 
time delays Tgk- Subsequently, the target location can be computed from the solution to a set of equations of the 
form (1), viz., 

Ttk = ^ ^\/ {xtk - xf + {ytk - yf + \J {xrt - xf + {yre. - yf^ ■ (8) 

The unknown complex ampUtude C is treated as a nuisance parameter in the estimation problem. 

Let the unknown target location X = {x,y) , unknown time delays delays rek, and unknown target complex 
ampUtude C = + iC^) where the notation specifies the real and imaginary components of (. 

We refer to the processing for estimating the target location as non-coherent or coherent. The received signal 
introduced in (7) is adequate for the coherent case, where the transmitting and receiving radars are assumed to 
be both time and phase synchronized. As such, the time delays information, r^^, embedded in the phase terms 
may be exploited in the estimation process by matching both amplitude and phase at the receiver end. In contrast, 
non-coherent processing estimates the time delays rik from variations in the envelope of the transmitted signals 
Sk it) . A common time reference is required for all the sensors in the system. In this case, the transmitting radars 
are not phase synchronized and therefore the received signal model is of the form: 

M 

n{t) = ^aekSk{t-Tek)+we{t), (9) 

fe=i 
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where the complex ampUtude terms a£k integrate the effect of the phase offsets between the transmitting and 
receiving sources and the target impact on the phase and amphtude of the transmitted signals. These elements are 
treated as unknown complex amplitudes, where a^k = o-fk ~^ J'^ik- define the following vector notations: 

a = [ail, Q!i2, aek, ckmn]'^, (10) 
= Re (a) ; =lm (a) , 
where Re (•) and Im (•) denote the real and imaginary parts of a complex- valued vector/matrix. 

III. Localization CRLB 

The CRLB provides a lower bound for the MSE of any unbiased estimator for an unknown parameter(s). Given 
a vector parameter 6, constituted of elements 9i, the unbiased estimate di satisfies the following inequality [19]: 

var(?,) > [J-1(^)].., i = l,2,... (11) 
where [J~^ (^)] j, diagonal elements of the Fisher Information matrix (FIM) J (6*). The FIM is given by: 



J (0) = Ee 



— logp(r|6') ( — logp(r|6') 



(12) 



where p {r\6) is the joint probabiUty density function (pdf) of r conditioned on 6. 
The CRLB is then defined: 

CcflLB = [JW]"'. (13) 

Sometime, it is easier to compute the FIM with respect to another vector ij:, and apply the chain rule to derive the 
original J [6) . In our case, since the received signals in both (7) and (9) are functions of the time delays, r^fc, and 
the complex ampUtudes, by the chain rule, J {6) can be expressed in the alternative form [19]: 

J(0)=PJ(V)P^, (14) 

where V is a vector of unknown parameters, and it incorporates the time delays. Matrix J (i/') is the FIM with 
respect to V, and matrix P is the Jacobian: 

From this point onward, we develop the CRLB for the case of non-coherent and coherent processing, separately. 

A. Non-coherent Processing CRLB 

For non-coherent Processing, there is no conmion phase reference among the sensors. Consequently, the complex- 
valued terms aik incorporate phase offsets among sensors and the effect of the target on the phase and complex 
ampUtude, following the definitions in (10). The vectors of unknown parameters is defined: 

9nc=[x,y,a'^,a'\' . (16) 
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The process of localization by non-coherent processing depends on time delay estimation of the signals observed 
at the receive sensors and also on the location of the sensors. To gain insight into how each of the factors affects 
the performance of localization, we utiUze the form of the FIM given in (14). We define the vector of unknown 
parameters: 

V'nc= [T,a^,a']^, (17) 

where a is given in (10) and r = [th, T12, r^fc, tmn^ ■ We are interested only in the estimation of x and y, 
while a^, act as nuisance parameters in the estimation problem. 

Given a set of known transmitted waveforms Sk{t — Tgk) parameterized by the unknown time delays r^j,, which 
in turn are a function of the unknown target location X = {x, y), the conditional, joint pdf of the observations at 
the receive sensors, given by (9), is then: 



1 ^ f 

P (r-li/Jnc) oc exp i / 

I O T 



M 



(t) - y^aekSk [t - nk) 



k=l 



dt 



The matrix P„c for (16) and (17), to be used in (14), is defined as: 



(18) 



p 



d _T 
dx ' 




_a_ 

dx 




d _T 


^ ("^)" 


d_ 

dy 




dr 






da' 


da"- 


9a" 




da" 


dr 


9a« 




da' 
da' 



{2MN+2)x3MN 

a dr 
where -^t is standard notation for taking the derivative with respect to x of each element of t, and — — =■ 

ox Q^H 

the Jacobian of the vector t with respect to the vector a^. The subscript denotes the matrix dimensions. 
It is not too difficult to show that using (8), the matrix P„c can be expressed in the form: 



(19) 



denotes 



(20) 



1 H2xMAf 02x2MAf 

^ O2MNXMN 'i-2MNy.2MN 

where is the all zero matrix, I is the identity matrix, and H e ji^'xmn incorporates the derivatives of the time 
delays in (8) with respect to the x and y parameters. These derivatives result in cosine and sine functions of the 
angles the transmitting and receiving radars create with respect to the target, incorporating information on the 
sensors and target locations as follows: 



H = 

The elements of H are given by: 



O-tXM "I" O-rxN 
bfXM "I" brxN 



(21) 
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(22) 



atx^ = cos 4>k; btxk = siri't'k; fc = l,..,M, 
ttrxe = cosLpi; brxe = sin (ft; i=l,..,N, 

(f>k = tan-i (y^iMh.) . = tan-i (m^) , 

^'^ \X-Xtk ) ' ^ \X — Xri ) ' 

where the phase is the bearing angle of the transmitting sensor k to the target measured with respect to the 
X axis; the phase (^^ is the bearing angle of the receiving radar i to the target measured with respect to the 
X axis. See illustration in Figure 1. For later use, we apply the following definitions: <^ = ^2, ^m]'^. 

= [</3i,</92,-,Viv] , ate = K btX2 ) •••) ^(Xm] 

and i>rx ~ [^rxi? ^rxg? *••) ^rxM] • 

An expression for the FTM J {tjjnc) , is derived in Appendix I, yielding: 



J(V'nc) = 



* nc 



(23) 



J {3MN)x(3MN) 

with the block matrices S„c, ^a, and V„c defined in the Appendix I in (92), (97), and (100), respectively. 
In order to determine the value of J {One) , we use (23) and (20) in (14), to obtain the following CRLB matrix: 



CcRLBnc = J ^ {One) = 



2/al 



HS„,H^ HV 

nc 

VLH^ a„ 



(24) 



The CRLB matrix is related to the sensor and target locations through the matrix H, and to the received waveforms 
correlation functions and its derivatives through the S„c and V„c matrices. 

1) Orthogonal Waveforms: When the waveforms are orthogonal, (92), (97), and (100) simplify to (101) in 
Appendix I. This simphfication enables to compute the CRLB (24) in closed form. We perform this calculation 
next. 



While the CRLB expresses the lower bound on the variance of the estimate of 9^ 



[x,y, 



a^, a^] , we are 



really interested only in the estimation of x and y. The amplitude terms and serve as nuisance parameters. 
For the variances of the estimates of x and y, it is sufficient to derive the 2x2 upper left submatrix [CcRLBnc\2x2 — 

"(J(^„c))"'" . 
L J 2x2 

Proposition 1: The CRLB submatrix [CcRLBnc\2x2 target locaUzation in the non-coherent case with orthog- 
onal signals is: 



[C 



Ci?LB„J2x2 ^ 



(HS„cH"^) 



Proof: From (101) in Appendix I, we have for terms of (24): 

S„e = 47r2/32 [diag(a)B diag(a*)] , 
V =0 

* nc ^1 

Aq, = 1.2MNX2MN- 



(25) 



(26) 
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In (26), diag(a) denotes a diagonal matrix with the elements of vector a. Matrix B = diag (l f^R,^])^ 
with denoting the normalized elements (Sn^ = /3ft//3, and 1 = [1, 1, ...1]^, 1 G R^^^. Using (26) in (24), it is 
easy to see that 



(,2 _i 

[CcijLs„J2x2 = TT^ (HS„cH'^) (27) 

9xnc ^nc 
f^nc 9ync 



c9ync - K 



nc 



where: 

n - .. 

M AT ^ 2 

k=li=l 

M w 2 (28) 

fe=l^=l 

/inc = - E E l"^fe| (<Jtefc + drxe) {hxk + ^rxj • 

fc=l£=l 

This concludes the proof of the proposition. ■ 
It follows that the lower bound on the variance for estimating the x coordinate of the target is given by 

9xnc9Vnc "'nc 

Similarly, for the y coordinate, 

^LcHB = Vnc ■ (30) 

The terms 9x„a' 9yna' ^nd hnc are sunnmations of atxk> drxe^ hx^ and brxe terms that represent sine and cosine 
expressions of the angles and ip, and therefore relate to the radars and target geometric layout. It is apparent 
that for the non-coherent case, the lower bounds on the variances (29) and (30) are inversely proportional to the 
averaged effective bandwidth /3^, and SNR = (see expression for ry„c in (28)). It is interesting to note that 
r]nc is actually the CRLB for range estimation in a single anterma radar, based on the one-way time delay between 
the radar and the target (see for example [19]). The other terms in (29) and (30) incorporate the effect of the sensors 
locations. 



B. Coherent Processing CRLB 

We recall that in the section on the signal model, we defined the complex amplitude aik associated with the path 
transmitter k target receiver £. In the non-coherent case, the complex amplitude is a nuisance parameter in 
estimating the target location x, y. In the coherent case, the transmitting and receiving radars are assumed to be 
phase synchronized. By ehminating the phase offsets, the signal model in (7) applies, and the nuisance parameter 
role is left to the complex target amplitude C = + iC^- The coherent approach to localization seeks to exploit 
the target location information embedded in the phase terms exp (— 27r/cT£fe) that depend on the delays r^fe, which 
in turn are function of the target coordinates x, y. 

Define the vector of unknown parameters: 
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(31) 



As before, define a second vector of unknown parameters in terms of the time delays t (rather then the target 
location), 

^-0= [t,C^C']^, (32) 

to be used in (14) to derive the CRLB. In comparing the coherent case in (32) with the non-coherent counterpart 
in (17), we note that V'nc incorporates the vectors and a^, while -^c is a function of the scalars C,^ and . The 
reduction in the number of unknown parameters is made possible through the measurement of the phase terms of 
and . 

For coherent observations, the conditional, joint pdf of the observations at the receive sensors, given by (7), is 
of the form: 



p(r|Vc) ocexpi -X-Y^ I 



M 



(t) - X^Cexp (-2-Kfcnk) Sk {t - Ttk) 



fe=i 



dt 



(33) 



We follow the same process used in Section III-A, to develop the CRLB for the coherent case based on the 
relation in (14). The matrix Pc takes the form: 



d0c 



H 0mNx2 
O2XMN I2x2 



(34) 



4x(MAf+2) 

where matrix H has the same form as in (21), since it is independent of the nuisance parameters in both cases. 
An expression for the FIM matrix, J (^c) > is derived in Appendix II, yielding: 



J(Vc) = 



(35) 



{MN+2)x{MN+2) 

where the submatrices are found in Appendix 11 as follows: Sc in (105), A^c in (108), and Vc in (111). 
The CRLB matrix for the coherent case is then found substituting (34) and (35) in (14) and (13), obtaining: 



CcRLBc 



(36) 



HScH^ HVc 
2M V^H^ A„, 

As in Section 111-A, we develop the closed form solution to the CRLB matrix in (36) for the case of orthogonal 
waveforms. Since we are interested only in the lower bound on the variances of the estimates of x and y, the 



submatrix [Ccrlb, 



is derived and evaluated next. 



Proposition 2: The CRLB 2x2 submatrix for the coherent case and orthogonal waveforms is: 

[CciiI,Bj2x2 =T772 



(HSeH^ - HVeA-^V^H^) 



(37) 
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Proof: From (112) in Appendix II we have the values of the matrices Sc, Aq,c, and Vc for orthogonal 
waveforms. Using this and H defined in (21) in (36), the CRLB matrix Ccrlb^^^ is obtained. Consequently, the 
submatrix [Ccrlb^]2x2 computed in Appendix III resulting in the form given in (37). 

This completes the proof of the proposition. ■ 

From (37) and (112), it can be shown that [Ccrlb^\2x2 '^^^ expressed as: 



[CckLbJ2x2 



Qxc he 
he 9vc 



(38) 



2 



where the various quantities are as follows: 

Stt^/KICIVt?,)' 
M N 2 / M N 

9xc = J2J2fRk {hxk + brxe) - WW S J2 (btxk + ^rxj 

fe=i^=i \k=ie=i 

M N 2 / M N \ ^ 

9yc = E E/«fc (otefc + cirxj - igjv E E + a™ J , (39) 

k=U=l \k=U=l / 

M N 

hc = - Y^HfRk iatxk + a-rxi) {hxk + brxi) 

k=U=l 
M N M N 

+ Am E E («fxfc + arxe ) E E ihxk + brxe ) ■ 
k=U=l k=U=l 

The lower bound on the error variance is provided by the diagonal elements of the [CciJiScor]2x2 submatrix 
and are of the form: 



The terms Qx^, gy^, and he are sunnmations of atxji. «rx^. hxk and brxe that represent sine and cosine expressions 
of the angles (p and ip and therefore relate to the radars and target geometric layout, multiplied by the ratio terms 
fuk = p')' Iiivokiiig tli^ narrowband signals assumption 01/ fe <C 1 it follows that fn^. ~ 1. These terms 
have some additional elements when compared with the non-coherent case. It is apparent that for the coherent 
case, the variances of the target location estimates in (40) are inverse proportional to the carrier frequency 

C. Discussion 
We make the following observations: 

• The lower bound on the variance in the non-coherent case is inversely proportional to the averaged effective 
bandwidth (3. For the coherent case, with narrowband signals, where /J^/Zc ^ 1. the locaUzation accuracy is 
inversely proportional to the carrier frequency fe and independent of the signal individual effective bandwidth, 
due to the use of the phase information across the different paths. It is apparent that coherent processing offers 
a target localization precision gain (i.e., reduction of the locaUzation root mean-square error) of the order of 



September 24, 2008 



DRAFT 



13 



/c//3, which we refer to as coherency gain. Designing the ratio /c//3 to be in the range 100-1000, leads to 
dramatic gains. 

• The term r]c in (39) is the range estimate based on one-way time delay with coherent observations for a radar 
with a single antenna [30]. 

• The CRLB terms are strongly reliant on the relative geographical spread of the radar systems vs. the target 
location. This dependency is incorporated in the terms Qx^^/x^i 9v^c/yc hnc/c- It is apparent from (40), 
(29) and (30) that there is a trade-off between the variances of the target location computed horizontally and 
vertically. A set of sensor locations that minimizes the horizontal error, may result in a high vertical error. For 
example, spreading the transmitting and receiving radars in an angular range of — (tt/IO) to +(7r/10) radians 
with respect to the target, will result in high horizontal error while providing low vertical error, as we would 
expect intuitively. This is caused by the fact that the terms Qxr^ahxc are summations of sine functions and 
dvrtchvc summation of cosine functions of the same set of angles. In order to truly determine the minimum 
achievable localization accuracy in both x and y axis, we need to minimize the over-all accuracy, defined as 
the total variance al = {<j1^crb + '^I^crb)- 

• The message of dramatic improvement in localization accuracy needs to be moderated with the observation that 
the CRLB is a bound of small errors. As such, it ignores effects that could lead to large errors. For example, 
MIMO radar with distributed sensors and coherent observations is subject to high sidelobes [IJ. Additionally, 
a phase coherent system is sensitive to phase errors. These topics are outside the scope of this paper, but they 
should be kept in perspective. 

• The lower bound as expressed by the CRLB, provides a tight bound at high SNR, while at low SNR, the 
CRLB is not tight. As stated in [33], the MLE is asymptotically unbiased and its error variance approaches the 
CRLB arbitrarily close for sufficient long observation time, with the condition that the MLE is not subject to 
ambiguilies. As the MLE of the time estimates is based on matched filters at the receiver end, the ambiguity 
features of the signal waveforms arise in low SNR conditions and predominate the estimation capabilities, 
causing erroneous time estimates. As the ambiguity problems are usually addressed trough the signal waveform 
design, a more rigid bound needs to be found for the localization variance in the low SNR case. 

IV. Effect of Sensors Locations 

The CRLB for target locaUzation with coherent MIMO radar shows a gain, i.e., reduction in the standard deviation 
of the locaUzation estimate, of fd P compared to non-coherent localization. Yet, the CRLB is strongly dependent 
on the locations of the transmitting and receiving sensors relative to the target location, through the terms gx„^/xc > 
9ync/vc ^nc/c- To gain a better understanding of these relations, and set a lower bound on the CRLB over all 
possible sensor placements, further analysis is developed in this section. 

We introduce the following general notation: for any given set of vectors ^ = (^i, ^2, --v^i) and {1^1, 1^2, kl): 
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i=l 
L 



(41) 



i=l 



The terms and in (28) can be expressed using the conventions defined in (41) and terms defined in 
Section III-B, viz.: 



MN 



T{hl)+T{hl,)-[T{h,x)f-[T{hr 



and 



9y^ = MN [t (e^,) + T (a^ J - [T (a,,)]^ - [T (a.,)]^] , 
where the narrowband signals assumption is applied. Similarly, the term he in (39) can be expressed: 



(42) 



(43) 



hc = MN[T{atxhtx)+T{airxhrr) (44) 
-T (atx) E (bte) - T {arx) E (b™)] . 

Since o^^^ +^Lfe — ^^^^ +sin^ 4>k — ^ ^rx< +^rx« — ^o^^ <^^ + sin^ (^^ = 1, the following conditions apply: 

T(afJ+T(bL)=l 



T(a2J+T(by =1 



0< [T(ai,)]' < 1 
0< [T(bt,)]' < 1 
< T (a2^) < 1 



< 1 



< [T (a^ 
0<T(a2^)<l 



(45) 



0<T(by <1; 0<T(by <1. 
We seek to find sets of angles (/)* and <^*, that yield sets of cosine and sine expressions a^^, a*^,, bj^, b*^ 
for which the values of the Cramer-Rao bounds for locaUzation along the x and y axes ((jI^crb ^'^'^ ^y^cKB' 
respectively) are jointly minimized, that is: 



nunimize 



(46) 



This is equivalent to minimizing the trace of the CRLB submatrix [CcrlBc\2x2- exphcit minimization problem 
is formulated introducing the objective function /q: 



minimize fo {^tx,B.rxMx,'^Tx) = glUl^-h'i 



(47) 



subject to constraints (45). 
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This representation of the problem is not a convex optimization problem. The next steps are undertaken in order 
to formulate a convex optimization problem equivalent to (47), i.e., a convex optimization problem that can be 
solved through routine techniques and from whose solution it is readily possible to find the solution to (47). 

In [28], it is shown that for a given positive definite matrix, in our case [CcrlbJ2x2' ^^'^ its inverse matrix F, 
in this case: 

F = l 

Vc 

the following relation exists between the diagonal elements of these matrices: 

[CcflLBj,,> » = 1,2. (49) 

EquaUty conditions apply for aU i iff F is a diagonal matrix, i.e., he = 0. Now, observe that the inverse of the 
elements on the diagonal of F are lower bounding the elements on the diagonal of the matrix CcrlBc for any 
atx, ^rx) btx, brx' We then define the objective function /o {atx, arx, btxj b^x) , and the optimization problem 



min /o (ate, arx, i>tx, Wx) = 




subject to (45). 

The new objective function and the original objective function are related as /o (a^^;, sl^x, btx, b^x) > /o (atx, arx,bta;, b^x), 
with equality for he = 0. Substitute the values of Qx^ and gy^ from (42) and (43) in the objective function of (50) 
to obtain 



dye 

-hr. 



-he 

Qxc 



(48) 



1/ {tIcMN) 

2 - T (bL) - T (b2 J - [T (a^)]^ - [T {^xt 

1/ {r]cMN) 



/o(ate,a,x,bte,b,x)- „ ™.,2 ^_^/K2 M2_rT^„ ^^^^ 



T (by + r (by - [T (bte)r - [T (b.xr 

It is apparent that the denominator of the first summand is bounded by: 

< 2 - r {hi) - T {h%) - [T {B.tx)f - [T Mf <2-T (hi) - T {h%) , (52) 
and the denominator of the second summand is bounded by: 



< T [hi] + T (b^x) - [T (btx)]' - [T {hrxf < T (hi) + T (hi,) . (53) 

'a convex optimization problem is of the form [32] 

minimize /o {x) 

subject to fi {x) < 
12 j O'j^j = 

for some constants Oj, i, j, i = l,...,m, j = l,...,p, and where /o,---,/m are convex functions. 
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Denote T (b^^) + T (b^^) = n, and let T (a^) = T (a^^) = T (b^) = T (b^^) = 0. Then, from (51)-(53) and 
(50), we obtain the following problem: 

— / N 1 1 

minimize /o (/x) = h - 

jii 2 — fj, fj, 

subject to M - 2 < (54) 

-H < 0. 

The objective function /o (/u) is still not convex. The epigraph form is a way to introduce a Unear (and hence 
convex) objective t, while the original objective /o is incorporated into a new constraint /o — i < 0. The key point 
here is that while /o is not convex, the constraint /o — t < can be transformed to a convex form. After some 
simple algebraic manipulations, the epigraph form turns into the following convex problem: 

minimize t 
subject to 

t/j,^ - 2i/i + 2 < ^^^^ 
yti- 2 < 
-M< 
-t<0 

A convenient way to solve this convex optimization problem is to employ the concept of Lagrange duality and 
exploit the sufficiency of the Kamsk-Kuhn-Tucker (KKT) conditions [32]. The Lagrangian of the problem in (55) 
is given by: 

t,\)=t + Xi {ti? - 2tiJL + 2) + A2 (/i - 2) - Asm - Mt, (56) 

where Aj, i= 1, .., 4 is the Lagrange multiplier associated with the ith inequality constraint /, [fi, t) < 0. 

The KKT conditions state that the optimal solution for the primal problem (minimization of f in (55)) is given 
by the solution to the set of equations: 

i, A) 



dL{iJL, t, A) 



= (57) 

= 



m 

fi{ti,t)<0; i = l,..,A 

Ai>0; j = l,..,4 

Xifi{ii,t) = 0; z=l,..,4. 
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(59) 



Applied to (55) and (56), these equations specialize to 

Ai (2i/x - 2t) + A2 - As = (58) 
1 + Ai (/i^ - - A4 = 
Ai (i/x^ - 2i/x + 2) = 

A2 (/I - 2) = 
-A3/X = 
-Xit = 0. 

It is not difficult to show that the solution to this system is given by 

t* =2 
A* = 1 

-^2 = ^3 — -^4 — 

RecalUng that ^ = T (bfj.) + T (b^^.) , the optimal solution can be rewritten as: 

li* =T{hll)+T{h*Z)=l. (60) 

In addition to (60), ci^^., a*^., b^^, b*^. have to satisfy the relations (45), and the equality conditions for (49), (52) 
and (53), viz., 

T(ar2) + r(a;2) = i 
T(b*j = 0; r(b;j = o 
r(a*j=:0; r(a;j = o 

Substituting these results in (42) and (43), we compute the optimal 5*^ and g*^ , 

gl = gl = MN. 

It follows that the minimum value of the trace of the Cramer Rao matrix [CcrlBc^^] , /o in (47), is given by: 



(61) 



The final step in determining the effect of sensor locations on the locaUzation CRLB is to recall that the 
multivariable argument of /o in (62) is actually a function of the transmitting sensors angles (f)i~, k = 1, . . . , M, 
and receiving sensors angles ipe, 1= 1,. . . ,N (see definitions in the previous section). What are then the optimal 
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sets (f)* and ip* that minimize the variance of the locaUzation error? The optimal angles can be found from the 
relations (61). For example, for the cosine of the transmitters bearings T (a^^) = 0, means 

1 

-^cos.^^ = 0. (63) 

A symmetrical set of angles of the form (p* = ^4>i\4>i = 4>„ + ^^^^g^!* = li --^M-.M 2|, is a solution to 
(63) for any arbitrary 0^. The same solution is obtained for the sines, T ih^^) = 0. The relations T(a*^) = 0, 
T (b*^) = lead to a solution constituted by a symmetrical set of angles (p* of the same form as </>*. The relation 
T (ajj-bj^.) + T {a.*^b*^) = expressed in terms of angles is 

1 ^ 1 ^ 

— COS sin (kk + j;^^ cos sin = 0. (64) 

fe=l £=1 

It can be shown that (64) is met by angles 0^ and (p^ symmetrically distributed around the unit circle, but the 
number of sensors has to meet M > 3, iV > 3. The condition T (b^J) + T (b*^) = 1 in (61), expressed in its 
explicit form, is 

1 ^ 1 " 

— cos' <l>l + j^Y. '^os' = 1- (65) 

The symmetrical set of angles that meet (63) and (64) provide jj cos' = jf cos' = | and there- 
fore meet the requirement of (65). The same applies to T (aj') -|-T (a*') = 1 , where we have Z^feli sin' (p^ = 
JV 2^^=ism iPi=^. 

We conclude that M > 3 transmitting, and N > 3 receiving sensors, symmetrically placed on a circle around 
the target at angular spacings of 2it/M and 2it/N, respectively, lead to the lowest value of the locaUzation CRLB. 

This result can be extended by noticing that relations (61) also hold for any superposition of synometrical sets 
containing no less than 3 transmitting and/or receiving sensors. Therefore, the complete set of optimal points is 
given by: 




;Z^>3;^Z„=M | 



(66) 



where the total number of transmitting (M) and receiving (iV) radars may be divided into V and U sets of 
symmetrically placed radars, each set consists of and Z^ radars, respectively. The angles 4>^ and are an 
initial arbitrary rotation of the symmetric sets Z^, and Z„, correspondingly. 

As a special case, it is interesting to evaluate the CRLB in (38) with 1 transmitter and MN receivers, i.e., a 
Single-Input Multiple-Output (SIMO) system. This scheme makes use of {MN + 1) radars instead of (M -|- N) 
radars used in a MIMO system with M transmitters and N receivers. From (66) it is apparent the this case does not 
provide optimaUty since the number of transmitters is smaller than 3. To evaluate ct^^crb + '^y^CRB for this setting 
we assume 1 transmitter is located at an arbitrary angle with respect to the target, and a set of MN receivers 
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are located symmetrically around the target, at angles y'* that follow the condition in (66). The expressions in (42), 
(43), and (44) reduce to the form: 



1 



(67) 



g,^=MN[T{hl,)-[T{hr.f 

gy^=MN[T{4,)-[T{Rr.)f 
he = MN [T (a,^b^^) - T (a^^) E (V^)] = 0, 

and the trace of the CRLB submatrix [CcrlBc]2x2- defined by /o (ata;,arx,btx,brx) = (^I^crb + ^leCRB = 

4??c 



3xc+3yc is 
9xc9yc-ht 



/o(at 



MN' 



(68) 



This result expresses an increase in the estimation error in the factor of 2 when compared with M transmitters and 
N receivers given in (62). 



A. Discussion 

The following comments are intended to provide further insight into the results obtained in this section. 

• From (62), the lowest CRLB for target localization utilizing phase information is given by 2ric/ (MN). We 
interpret the reduction of the CRLB by the factor MN/2 compared to a single anteima range estimation given 
by r]c as a MIMO radar gain. This gain reflects two effects: (1) the gain due to the system footprint; (2) the 
advantage of using M transmitters and N receivers, rather than, for example, 1 transmitter and MN receivers. 
The latter gain is apparent when MN :$> {M + N). 

• The CRLB obtained through the use of a single transmit anteima and MN receive antermas in (68) is 
477c/ (MN). It follows that MIMO radar, with a total of M + N sensors, has twice the performance (from the 
point of view of localization CRLB) of a system with a single transmit anteima and MN receive antennas. 

• The best accuracy is obtained when the transmitting and receiving radars are located on a virtual circle, centered 
at the target position, with uniform angular spacings of 2tt/M and 2'k/N, respectively, or any superposition 
of such sets. 

• The optimization analysis presented in this section is intended to provide insight into the effect the sensors 
locations have on the CRLB. Naturally, in practice, it is not possible to control in real time the location of the 
sensors relative to a target. However, the results here teach us that selecting among the sensors those who are 
most symmetrical with respect to the target may lead to the most accurate locaUzation. 

So far we have focused on the theoretical lower bound of the localization error. In the next section, we discuss 
specific techniques for target localization and their performance as a function of sensors locations. For this purpose, 
the GDOP metric and GDOP contour mapping tools are introduced. 
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V. Methods for Target Localization 

In Section III, was formulated the lower bound on the variance of any localization estimate. Here, it is of interest 
to discuss some specific target localization estimators. In particular, two estimators are presented: the MLE and the 
BLUE. The MLE is motivated by its asymptotic optimahty, while the BLUE by its closed form expression. 

A. MLE Target Localization 

The MLE is a practical estimator in the sense that its application to a problem of observations in white Gaussian 
noise is relatively straightforward. Moreover, under mild conditions on the probability density function of the 
observations, the MLE of the unknown parameters is asymptotically unbiased, and it asymptotically attains the 
CRLB [19]. 

For the case of coherent MIMO radar, the signal waveform received by radar £ is given in (2). The MLE of the 
unknown parameter vector 9 = [x,y,(\^ given the observation vector r is given by [19]; 

^ML = arg \ max [logp {r\0)] \ , (69) 



where p {r\9) is given by (33) noting that the time delays rek are known functions of x and y. To jointly maximize 
logp (r|^) with respect to ^ = [x, y, C,]^ , we start by maximizing it with respect to (: 

^logp(r|x,y,C) |^=f=0. (70) 

Using (33) in (70), the estimate C can be found, and it is a function of x and y. By substituting it back into (69), 
it is said to compress the log-Iikelihood function [31] to logp (t\x, y, . The MLE of the target location is then 
given by 

■^logp(r\x,y,C^ \^^Sml=0 

^logp(r|a;,y,c) |j/=?Mi.= 0. (71) 

Since a closed form expression can not be found for the MLE in (71), numerical methods need to be applied. A grid 
search or an iterative maximization of the likelihood function needs to be performed to determine xml and yML- 
This might involve a significant computational effort. In practice, we can limit the search grid for high resolution 
target localization estimation to an area around a coarse initial estimate obtained by the non-coherent approach. 

B. BLUE Target Localization 

The MLE presented in Section (V-A) does not lend itself to a closed form expression, and numerical methods 
need to be used to solve it. A closed form solution to the target localization can be obtained by application of the 
BLUE. 

To formulate the BLUE, it is necessary to have an observation model in which observations change linearly with 
the target location coordinates. That is because it is inherent to the BLUE that the estimate is linear. To this end. 
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we formulate a model in which the time delays are "observable." Let the observed time delay associated with a 
transmitter-receiver pair be /i^fc, then 

m = rek + £ek, V/c = 1, .., M, Z = 1, .., iV, (72) 

where Sik is the "observation noise." In practice, the time delays are not directly observable. Rather, they are 

estimated, for example by maximum likelihood, from the received signals. Then, the term eek is the time delay 

estimation error. Our BLUE estimation problem of the target location should not be confused with the estimation 

of the time delays. The estimation of the time delays is just a preparatory step in setting up the "observations" of 

the BLUE model. Once, the observation model has been set up, it is necessary to ensure that the model between 

the time delays and target location is linear. Setting the origin of the coordinate system at some nominal estimate 

of the target location, and preserving only Unear terms of the Taylor expansion of expressions such as in (1), we 

can express the time delays as Unear functions of x and y, 

X y 
Tek ^ - - (cos (pk + cos (fi) - - {smcpk + sm(pe) , (73) 

where the angles and ipt are the bearings that the transmitting sensor k and receiving sensor £, respectively, 
subtend with the reference axis (with the origin at the nominal estimate of the target location). Note that the 
definitions of the angles here are a little different than the angles defined in Section III and also denoted <p and 
if. Here, the vertex of the angles is an arbitrary point in the neighborhood of the true target location. In Section 
in, the vertex is at the true target location. Since only the vertex is different, we preserved the same notation for 
simpUcity sake. UtiUzing definitions (22), we can express the linear model in the following simpUfied form: 

nk = -- {atxk + O-rxe) - - {hxk + ^rxj • (74) 

Letting, r = [ni, ri2, tmn]^ and the vector of unknowns 6 = [x, y, C]^, we write (74) in vector notation as 
follows: 



T = m, 

where the angle dependent matrix D is defined as: 



D 



The observation model (72) can then be expressed as 

At = D6I + e, 



(75) 



(76) 



MJVxS 



(77) 



where /i — [/in, /ii2, i^mn]^ , and e = [en, £12, £mn\^ is the MN x 1 observation noise vector. To reiterate, 



a key difference between the MLE and BLUE models is that the MLE target locahzation is carried out utilizing 
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signal observations (which are not linear in a;, y), while according to (77), the BLUE's "observations" are in the 
form of time delays. So an intermediate step of time delay estimation is implied. The time delays estimates used 
as observations can be derived for example by MLE as follows: 



/i^fe = arg max 



exp {j2'Kfcv) / re (i) si {t - v) dt 



(78) 



where v is a dunmiy variable for the time delay. 

We still need some characterization of the "noise" terms sek- It is shown in Appendix IV, that the maximum 
likelihood time delay estimates are unbiased with error covariance matrix 

1 



l-MNxMN, 



(79) 



where previous definitions of the various quantities apply. For the Unear and Gaussian model in (77), the BLUE 
is computed from the Gauss-Markov theorem [19] that states the BLUE of the unknown vector 6 is given by the 
expression: 



?5 = (D^C,-iD)-^D^CrV 
The theorem also establishes that the error covariance matrix is 



(80) 



Cs = (D^Cr^D)-\ 



(81) 



Using the time error covariance matrix and the linear transformation matrix D in (76), the following estimate 
for the target locahzation is obtained: 



X 

y 



-cGf 



M N 

EE(«t ) M^fc 

k=lf=i 
M N 

k=U=l 



(82) 



where fx^k are the time observations, and the matrix Gs is of the form: 



9ib92b - hi 



9 IB 

^B 921 



(83) 



The elements of matrix G^ are: 
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M N , / M N ^ ^ 



= EE (^*-^ + ^r.f - EE + ^™^) ) ' (84) 

k=ie=i \fe=i£=i 

M AT ^ / M N 

92B = EE + "™^)^ - MN ^ 
k=it=i \k=ie=i 



M N 



hs = -EE + '^™^) + ^™«)) ' 

fc=l^=l 

^ M N M N 

+ MivEE + "™^) EE + ^™^) ■ 



/c=l£=l fc=l^=l 

Using these results in (81) provides the MSE for the BLUE as follows: 



for the estimation of the x coordinate, and 



(85) 



,2 



for the estimation of the y coordinate. 

C. Discussion 

The following points are worth noting: 

• The BLUE estimator in (80) and its variance in (85) and (86) are provided in closed form. This enables analysis 
without extensive numerical computations. 

• In general, the variances (85) and (86) have similar functional dependencies on the carrier frequency and on 
the sensor deployment as the CRLB (40). The terms atx^, drxe^ hxk ^rx^ embedded in (85) and (86) 
relate the sensors layout to the variance of the BLUE . 

From the expressions of the variance of the BLUE, one can not readily visuaUze the effect of the sensors layout. 
A mapping method, acting as a design and decision making tool for MIMO radar systems, is proposed and evaluated 
in the next subsection. 

D. GDOP 

In Section IV, we discussed optimal sensor location for minimizing the CRLB. In practice, we are faced with a 
specific deployment of sensors, and we ask what is the localization accuracy for a given location of the target. GDOP 
is a metric that addresses this question. The GDOP is conraionly used in GPS systems for mapping the attainable 
locaUzation accuracy for a given layout of GPS satellites positions [28], [29]. The GDOP metric emphasizes the 
effect of sensors locations by normaUzing the localization error with the term contributed by the range estimate. 

The GDOP metric for the two dimensional case is defined: 
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GDOP = , (87) 

where cr^ and are the variances of locahzation on the x and y axis, respectively, and ct^ is the standard deviation 
of the time delay estimation error, assumed the same for all sensors. Inherently, the GDOP provides a normalized 
value that measures the relative contribution of the radars' location to the overall accuracy. When the BLUE is 
used, and the hnearity conditions hold, cr^ and are given by (85) and (86), respectively. Using the result in (79), 
CCTe for the time delay variance, we get the following GDOP expression: 



GDOPb = J^i^^^. (88) 

y .9iB.92B - K 

The GDOP reduces the combined effect of the locations of the sensors to a single metric. Once we get the values 
mapped, the actual localization error is easily derived by multiplying the GDOP value with ca^. 

Figure 2 and 3 present contour plots of the GDOP values for 3x4 and 7x7 MIMO radar systems, re- 
spectively. The sensors are positioned symmetrically around the origin. In Figure 2, the transmitting sensors 



are located at bearings = 

V = 



_ 27r(i-l) 



7 = 1 

3 ' ' 



, and the receiving sensors are positioned at bearings 



In Figure 3, the M = 7 transmitting sensors are positioned as a superposition 
of two symmetrical constellations: the first set includes three radars and the second four. The sets are located 
at bearings (j) = 0i = jg + i ~ 1, 3; 0i = f + ^^^^3—^, * = 4, 7 . The receiving radars, for 

this case, are set in a single symmetrical constellation with bearings if — (pi — ^ i = 1, 7 . The first 

noticeable factor in the comparison of the two plots is the higher accuracy obtained with seven radars compared 
to four radars. For example, the lowest GDOP value in Figure 2, for the 3x4 system is 0.4082, while with seven 
radars (see Figure 3), the lowest GDOP is 0.2020, corresponding to a 50% reduction. When a target is located 
inside the virtual {N + M)-sided system footprint, a higher localization accuracy is obtained than when a target is 
outside the footprint of the system. In particular, the best localization is obtained for a target at the center of the 
system. The increase in GDOP values from the center to the footprint boundaries is slow. Outside the footprint, 
the GDOP values increase rather rapidly. 

In Figure 4 and Figure 5, contours of seven non- symmetrically positioned radars are drawn. When the radars 
are relatively widely spread, as in Figure 4, there are stiU some areas with good measurement accuracy, though the 
coverage is shrunk compared to the case with symmetrical deployment of sensors in Figure 3. When the viewing 
angle of the target is very restricted, as in Figure 5, there is a marked degradation of GDOP values. 

These examples demonstrate the main theoretical result of Section IV, namely that a symmetrical deployment of 
sensors around the target yields the lowest GDOP values. Furthermore, calculating the lowest attainable GDOP value 



using the optimal results in (62) for a M x A?^ MIMO radar, we obtain a GDOP value of ^J^jMN, and for M = AT 
it is equal to ^/2jN^. As a numerical example, the lowest GDOPs in Figures 2 and 3 are -\/2/3 • 4 ~ 0.4082 and 



\/2/T^ ~ 0.2020, respectively. Comparing this with the results obtained in [29J for the case of passive GPS based 
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systems, with N satellites optimally positioned around the target, for which the lowest achievable GDOP value is 
2/\/]V, the MIMO system advantage is clearly manifested. 



In this paper, we have developed analytical expressions for the estimation errors of coherent and non-coherent 
MIMO radar using the CRLB. It was shown that when the processing is coherent and the phase is processed, 
there is a reduction in the CRLB values (standard deviation of the estimates) by a factor of /c//3 over the case 
when the observations are non-coherent. We referred to this gain as coherency gain. Expressions for the CRLB 
capture also the impact of the sensors geometry. Further minimization of the locaUzation error reveals a MIMO 
radar gain directly proportional to the product of the number of transmitting and receiving radars. The smallest 
CRLB is achieved when the transmitting and receiving sensors are arrayed synometrically around the target or any 
a superposition of such sets. The GDOP metric and mapping were introduced as a general tool for the analysis of 
the localization accuracy with respect to the given radars and target locations. These plots could serve as a tool 
for choosing favorable radar locations to cover a given target area. While localization by coherent MIMO radar 
provides significantly better performance than non-coherent processing, it faces the challenge of multisite systems 
phase synchronizing, and needs to deal with the ambiguities stenaming from the large separation between sensors. 

Appendix I 
Derivation of the FIM in (22) 

In this appendix, we develop the FIM for the unknown parameter vector ipnc, based on the conditional pdf in 



VI. Conclusions 



(18). The expression for J (tp) = E V^, logp {r\il)) (V^ logp (r|V')) = -E 



logp(r|^) 



is derived using: 




(89) 



e,e' = 1,..,N; k,k' = 1,..,M; 



The first derivative of p {r\ipnc) with respect to the elements of r is: 




(90) 



fe'=i 



d[sk{l 



dTik 



Ilk 



1 



dt. 
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Applying the second derivative to (90), define a matrix S„c with the following elements: 

^2 



[S. 



= E 



— / [aekSk{t-nk)(x^k'sl>{t-Tek') 

t'k' J 



Using matrix notation for compactness, 



drekdTi 

+ o^ek^k (t - Tik) aik'Sk' (t - Tik')] dt} 
92 



dnkdn 



— J Sk{t-nk)sl,{t-Tik')dt I 



S„c = -Q-2 Re [diag(a)Rs diag (a*)] , 
where diag(-) denotes a diagonal matrix, a was defined in (10), and we abuse the notation and let 



d 



,., dTikdri'k' 



The elements of matrix are defined as: 



S Skit- Tek) si, {t - Tek' ) dt 



£^e' 
The second and third terms in (89) define a matrix with the following elements: 

[Aa]ii' = [^a]{MN+i),{MN+i') = [J {i'nc)](MN+i),(MN+i') = ('^^"c)](2MJV+i),(2MAr+i') 



X^Sfe - rek) a}k'S*k' [t - nk') 



M 



.fc' = l 



M 



■'^sl {t - Uk) cx-tk- Sk- [t-Tek') 



fe'=i 



dt 



Re{[R.],,,}, 



and 



[Aa;]i,(MiV+i') = [Aa;](MJV+i),i' = [J (V'Tic)](MAr+i),(2MJV+i') = ["^ {i'nc)](^2MN+i),{MN+i') 



M 



d 



M 



(-j) (i - a^fe'Sft' - Tek') 



k'=l 



dt 



= -Im{[R,]..,}. 



In matrix notation. 



A„ = 



Re [Rs] - Im [Rs] 
- Im [Rg] Re [R^] 
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The fourth and fifth terms in (89) define the matrix V„c with the following elements: 



and 



d d 



: E 



— / [oiekSk{t-Tek)oilysl,{t-nk') 

t'k' J 



+oiik'^k (t - Tek) atk'Sk' {t - Tck')] dt} 
d 



= Re <^ aik 



(7^ C7^ 
[V„c]i,(MJV+j') = [J {i'nc)](^2MN+i),i' = [J {i'nc)]i^(^2MN+i') 



2 

: E 



d 



dTfk da 



d f 

—r — / [o:ekSk{t-nk)a*ek'S*k'{t-Tek') 

^ i'k' J 



+a|feSfc (t - Tek) aik'Sk' {t - Uk')] dt) 
d 



= - Im <^ 



'dTik 



In matrix notation: 



^ Re [diag(a)Rs] ; --^ Im [diag(a)Rs 



(98) 



(99) 



(100) 



Orthogonal Waveforms 

Orthogonality implies that all cross elements J Sk {t — r^k) s^, (t — Tuk') dt = 0, for I' and k ^ A;', and after 
some algebra, the matrices defined by (91)-(99) take the following form: 



[S. 



47r2/32 



[Aa]ii' — [A„](MAr+i),(MAr+i') - 
nc\ii' ~ 




(101) 



[V, 



0. 



Appendix II 
Derivation of the FIM in (34) 

In this appendix, we develop the FIM for the unknown parameter vector i^j^, based on the conditional pdf in 



(33). The expression for i{i)) = E | logp (r|V') (V^, logp (r|V'))^} = -E 



d logp(r|i/i) 



is derived using: 
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[Hi>c)],,=-E 



9 logp(r[V>e) 



[J(V'c)] 



(MJV+1),(MJV+1) 



c^J(MJV+2),(MJV+2) 



-E 
-E 



d logp(r|t/;c) 
logp(r|-0e) 



(MN+2),(MN+1) ~ E 

log p(r\tpc) 



[J (V'c)](MJV+l),(MAr+2) — [J (''/'c)] 
[J (V'c)]i,(MJV+l) = [J (V'c)](MJV+l),i' = 
[J (V'c)]i,(MJV+2) = [J (V'c)](MJV+2),i' = -E 

i = {£-l)*M + k, i' = {£' -l)*M + k', 
£,£' = 1,..,N; k,k' = l,..,M. 

The first derivative of p {r\ipc) with respect to the elements of r is: 



logp(r|V)e) 



3 logp(r|V>e) 



(102) 



d [logp(rlV'c)] _J_ r 

„2 J 



M 



(i) - ^ C exp {-j27rfcnk' ) Sfc' - rek' ) 



c 



H. 19 [ exp {j2wfcTik) si {t - Tek)] 
drek 



+ 



M 



fe'=i 



9 [exp i-j2Tr fcTik) Sk {t - nk) 
dTik 



Applying the second derivative to (103) define a matrix S„c with the following elements: 



[Sc 



(103) 



dt. 



(104) 



E I / [CC* exp (j27r/c (r^fe - r^'fe')) Sfe' - T^fe') Sfc {t - Tek) 

I OTlkOTlik' J 



CCexp (-j27r {rik - nk')) 5%, {t - Tek') Sk {t - Tek)] dt} 
52 



= Re Id 



dTekdTe'k' 



(exp {-j2-Kfc {Tek - Te'k')) [^slw) 



In matrix form, 



o2 

= ICI' Re {diag(e)R, diag (e*) } , (105) 

where the operator ^ and the matrix Rg were defined in Appendix I, e = [exp (— 27r/crii) , exp (— 27r/cTi2) , 
...,exp (-27r/cTAfAr)]. 
The second and third terms in (102) define a matrix Aac with the following elements: 



acjll 



= [A, 



E 



ac\22 — [J (V'c)](MJV+l),(MJV+l) 



[Hi'c)i 



N M 



2 I- \Ti./j(MJV+2),(MJV+2) 

M 

^ exp {-j2wfc {Tek - Tek')) Sk {t - Tek) si, {t - Tek') 



(106) 



EE/ 

M 

+ ^ exp (j27r/c {Tek - Tek')) si {t - Tek) Sk' {t - Tek') 
k'=i 

{N N M M ~j 
^^^^exp(-j27r/e {Tek Te'k')) [R-slii' ( ; 
e=U'=lk=lk'=l J 



September 24, 2008 



DRAFT 



29 




September 24, 2008 



DRAFT 



30 



Orthogonality implies that all cross elements / Sk{t — Tgk) s^, {t — Te'k')dt = oior £' and k ^ k'. Therefore, 
the matrices defined by (104)-(110) take the following form: 

i^i' 



[ i^i' (112) 
[Aac„j21 = [Aa„Jl2 = 
[V,„J,i = 27rCVc 

[Vc„Ji2 = -27rC^/c. 

where /ij^ = ^1 + When we invoke the narrowband assumption 0^/ fc ^ 1 it follows that /^^ ~ 1. 

Appendix III 
Computation of (36) 

The submatrix ['OcrlBc\2x2 defined as: 

[CciJLBe]2x2 = [J (^c)]2x2 • (113) 

For a given matrix of the form: 



HScH^ HVc 

where A^c is a diagonal matrix of the form A„c = d^2x2, and is some constant. 



(114) 



By definition, the value of [J (^c)]i i is obtained by: 



[J(^c)]7l 



J(^c) 



ea;(l,l) 



|J(ec 



(115) 



where |-| denotes the determinant, and J (^)ex(i,i) i^ ^ submatrix, obtained by removing the first row and the first 
colunm of the J {6c) matrix. The determinant of J (^c). using the property that the determinant of a matrix does 
not change under linear operations, is: 



|J(^c)| = 

This can be calculated and expressed as: 



HS.H^ - Vj^H^Ar^HV, 



(116) 



|J(^c)| = 



HScH^ - VTH^AJ^HV 



(117) 



Repeating the same for the matrix J (^c)ea;(i,i)- 
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J (^c)ex(l,l) = 

Using the same matrix manipulation, we get: 



HScH'^gj.n ,1) HV, 



(1-1) 

ex(,l) 



cex{l,) 



(118) 



HScH^ - V^H^A-iHVc 



J(^c)e.(i,i) = HSeH^ - V^H^A-iHV, 
and using terms (117) and (119) in (115) yields: 

[J(^c)]r,} = 

By definition, this expression is identical to: 

[J {e,)]-^ = [(HS,H^ - V^H^A-iHVe)"' 
Repeating the process for term located at (1,2), (2, 1), and (2,2), results in: 



1,1 



(119) 



(120) 



(121) 



(122) 



Appendix IV 

Derivation of Covariance of Observation Noise (78) 

For a set of received waveforms {t) , 1 < £ < N, the time delay estimates fi = [/Uii,/ii2, ■■■,IJ'Mn]^ are 
determined by maximizing the following statistic: 



liik = arg max 



exp(j27r/cw) / re (t) si {t - v) dt 



Equivalently, 



dv 



exp {j2Trfcv) j n {t) si {t - v) 



dt 



= 0. 



(123) 



(124) 



The time delay estimates are expressed in (72). The properties of the noise e^k can be computed from (8), and (2). 
It is not difficult to show that the following relation holds: 

dg{v) 



dv 



■ nek = 0, 



where 



and 



9{v) = C j exp [j2iTfc {v - Tik)] Sk {t - Tek) si {t - v) dt, 



nek = I —we{t)sl {t - v) exp {j2-Kfcv) dt. 



(125) 



(126) 



(127) 
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We wish to write (125) in the form of (72). With a few algebraic manipulations, including expanding g{v) in a 
Taylor series around rck, and neglecting terms o {rek —Tik f , it can be shown that 

m = Tek-\ 7 • (128) 

Air^f? (l + 1) C 

Comparing this with (72), and invoking the narrowband assumption Pl/f^ ^ 1, we have for the error term 

^ nek 

To find the first and second order statistics of e^k, we need the statistical characterization of nek- As previously 
stated, we assume the receiver noise we {t) is a Gaussian random process with zero mean and autocorrelation function 
a^S{T). Since nek is a linear transformation of the process we{t), since the mean we{t) is zero, E[nek] = 0. 
Similarly, it can be shown that 

, _ y£k ^ nm 

E[nekKJ = { „ „ .„ • (130) 

y£k = nm 




Using these results, we finally get 

leTT^iciVc" 



E [Ukf-nml = , ,^,2 .4 (131) 
/ c 

V^fc ^ nm 

concluding that the covariance matrix of the terms eek is given by: 

= 7j ^MNxMN- (132) 
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Fig. 1. MIMO radar system layout. 
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X axis 

Fig. 2. GDOP contours for a symmetric positioning of radai's around the axis origin: case (a) witli M=3 transmitting radars in a symmetric 
constellation with the receiving radars set organized in a symmetric constellation of N=4. 




X axis 

Fig. 3. GDOP contours for a symmetric positioning of radars around the axis origin: case (b) with M=7 transmitting radars in two symmetric 
constellations, one of 3 radars (in red) and second of 4 radars (in green). The receiving radars set are organized in a symmetric constellation 
of N=7. 
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X axis 



Fig. 4. GDOP contours for an asymmetric constellation of the radar set with M=7 transmitting radars and N=7 receiving radars. 




Fig. 5. GDOP contours for an asymmetric constellation of the radar set with M=7 transmitting radars and N=7 receiving radars, in the case 
where the radar are almost aligned. 
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